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Abstract 

A  three-dimensional,  steady-state  mathematical  model  is  described  to  investigate  the  fluid  flow,  species  transport  and  electrochemical 
reaction  in  the  PEM  fuel  cells  with  conventional  and  interdigitated  flow  fields.  The  multidimensional  characteristics  of  flow,  species  and 
current  distributions  are  computed  by  method  based  on  volume-control  finite-discrete  technique.  After  comparing  the  absolute  value  of 
convective  and  diffusion  component  of  oxygen  flux  quantitatively,  we  found  that  forced  convection  transport  mechanism  is  dominant  for 
the  interdigitated  flow  field  design,  however  the  diffusion  transport  is  dominant  for  the  conventional  flow  field  design.  The  performance 
and  the  pressure  loss  of  these  two  different  designs  are  calculated  and  compared;  results  show  that  interdigitated  flow  field  design  have 
better  performance  for  its  advantage  in  mass -transport  ability,  however  have  larger  pressure  loss  for  its  flow  through  the  electrode.  Finally, 
experimental  results  reported  in  the  literature  and  predicted  polarization  curves  are  compared  to  evaluate  the  numerical  model  employed. 
©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC) 
engines,  using  polymer  membrane  as  electrolyte,  are  con¬ 
sidered  to  be  promising  power  sources,  especially  for 
transportation  applications.  This  fuel  cell  has  many  im¬ 
portant  attributes  such  as  high  efficiency,  clean,  quiet,  low 
temperature  operation,  capable  of  quick  start-up,  no  liq¬ 
uid  electrolyte  and  simple  cell  design.  However,  before 
this  system  becomes  competitive  with  the  traditional  com¬ 
bustion  engine,  its  performance  and  costs  must  be  further 
optimized.  In  order  to  enhance  its  performance,  a  new  flow 
field  design  called  the  interdigitated  flow  field  has  been 
developed  and  proven  to  be  very  effective  by  experiment 
[1-3],  as  schematically  shown  in  Fig.  lb.  Comparing  to 
flowing  over  the  surface  of  the  electrodes  in  the  conven¬ 
tional  parallel-channel  flow  field  shown  in  Fig.  la,  these 
interdigitated  flow  fields  in  effect  have  converted  the  trans¬ 
port  of  the  reactant/product  gases  to/from  the  catalyst  layers 
from  a  diffusion  dominated  mechanism  to  a  forced  convec¬ 
tion  dominated  mechanism  by  forcing  the  reactant  gases 


*  Corresponding  author.  Tel.:  +86  571  8795  1764; 
fax:  +86  571  8799  1863. 

E-mail  address:  fanjr@mail.hz.zj.cn  (J.  Fan). 

0378-7753/$  -  see  front  matter  ©  2004  Elsevier  B.V.  All  rights  reserved. 
doi:10.1016/j.jpowsour.2004.05.010 


to  flow  into  the  electrode  in  order  to  exit.  The  transport 
rate  of  the  reactive  gases  to  the  catalytic  reaction  layer  is 
enhanced  and  the  electrode-flooding  problem  in  the  cath¬ 
ode  is  reduced,  thus  increasing  the  current  density  at  a 
given  cell  potential  and  the  limiting  current  density  signif¬ 
icantly. 

In  the  last  decade  a  number  of  one-  and  two-dimensional 
mathematical  models  [4,5]  have  been  presented  to  investi¬ 
gate  PEM  fuel  cells  with  conventional  flow  field  and  their 
performance.  These  models  cannot  account  for  effect  of 
the  flow  channel  width  on  the  velocity  distributions  and 
the  species  distributions.  More  recently,  three-dimensional 
models  [6,7]  based  on  the  computational  fluid  dynamics  ap¬ 
proach  have  been  presented  and  computed  by  commercial 
software  such  as  FFUENT  and  CFX.  However,  for  the  inter¬ 
digitated  flow  field  design,  some  studies  focusing  on  devel¬ 
oping  two-dimensional  mathematical  model  [8,9]  have  been 
conducted  to  delineate  the  flow  and  transport  mechanism. 
These  two-dimensional  models  cannot  consider  the  compo¬ 
sition  changes  along  the  flow  path,  thus  their  applicability  is 
limited.  Um  et  al.  [10]  carried  out  a  three-dimensional  com¬ 
putational  study  for  PEM  fuel  cells  with  conventional  and 
interdigitated  flow  fields.  Fully  three-dimensional  analysis 
of  transport  and  reaction  in  PEM  fuel  cells  with  the  inter¬ 
digitated  flow  field  was  presented. 
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Nomenclature 

C  molar  concentration,  mol  cm-3 

Cf  fixed  charge  concentration,  mol  cm-3 

D  diffusion  coefficient,  cm2  s_1 

F  Faraday  constant  =  96487  C/mol 

H  cell  height,  cm 

I  mean  current  density,  A  cm-2 

J  transfer  current  density,  A  cm-3 

kv  hydraulic  permeability,  cm2 

k§  electrokinetic  permeability,  cm2 

L  catalyst  layer  width,  cm 

N  component  number  and  electron  number 

participating  in  the  electrochemical  reaction 
P  pressure,  atm 

R  gas  constant  and  production/consumption 

rate  of  material 

S  stoichiometric  coefficient 

T  cell  temperature,  K 

u  velocity  vector,  cms-1 

Zf  fixed  site  charge 

Greek  letters 

a  transfer  coefficient 

y  concentration  parameter 

s  porosity  of  the  porous  media  material 
p  fluid  viscosity,  kg  cm-1  s-1 

p  density  of  the  fluid,  mol  cm-3 

a  conductivity,  £2_1  cm-1 

(p  potential,  V 

Superscripts 

eff  effective  value 

Subscripts 
a  anode 

c  cathode 

f  membrane  species  (the  fixed  charge  site  species) 
k  component 

m  membrane 

ref  reference 


The  major  objective  of  this  work  is  to  develop  a  three- 
dimensional  model  for  the  simulation  of  PEMFCs,  which 
can  be  applied  to  investigate  the  multidimensional  distri¬ 
butions  of  current  densities,  reactants  concentrations  in 
the  fuel  cell  with  the  interdigitated  and  conventional  flow 
fields,  knowledge  the  fuel  cell  work  principles  and  mass- 
transport  schemes  of  these  two  cell  design  in  more  detail. 
At  the  end,  polarization  characteristics  for  conventional 
and  interdigitated  flow  fields  are  predicted  and  compared. 
Estimated  and  experimental  results  in  the  published  lit¬ 
erature  were  compared  to  verify  the  mathematical  model 
used  in  this  paper,  reasonable  agreement  was  achieved.  In 


this  first  step,  water  transport  in  the  membrane  electrode 
assembly  is  not  modeled  in  detail.  The  effect  of  relative 
humidity  of  reactant  gases  on  the  performance,  and  the 
two-phase  phenomena  under  high  current  density,  is  not 
discussed  in  this  paper,  which  will  be  considered  in  the 
future  work. 


2.  Mathematical  model 

The  scheme  of  a  PEMFC  is  shown  in  Fig.  1.  Fig.  la  and  b 
is  flow  path  structure  for  the  conventional  and  interdigitated 
flow  fields,  respectively.  For  the  conventional  flow  field, 
two  flow  channels  around  the  current  collector  are  same, 
the  gas  mixture  flows  in  from  the  inlet,  and  flows  over  the 
electrode  surfaces  in  parallel,  and  then  out  from  the  outlet. 
For  the  interdigitated  flow  field,  however,  two  flow  channels 
around  the  current  collector  have  a  dead-end,  which  force 
the  gas  mixture  flow  into  the  backing  layers,  and  then  flow 
out  from  the  outlet  of  the  outlet  channel.  The  common  PEM 
fuel  cell  can  be  divided  into  seven  segments  across  the  fuel 
cell,  as  illustrated  in  Fig.  lc.  The  two  outer  segments  are 
graphite-plate  current  collectors  with  flow  channels  for  re¬ 
actants  and  products  to  enter  and  exit  the  cell,  respectively. 
The  adjacent  segments  are  electronically  conducting  porous 
backing  layers  that  allow  for  even  distribution  of  reactants 
to  the  anode  and  cathode.  The  centric  polymer  membrane 
spatially  separates  the  fuel  cell.  The  catalyst  layers  between 
the  membrane  and  the  backing  layers  are  domain  for  elec¬ 
trochemical  reactions.  The  catalyst  layers  provide  electron, 
proton  and  reactants  channels.  In  this  work,  the  computa¬ 
tional  domain  consists  of  three-dimensional  complete  fuel 
cell. 

The  assumptions  used  in  developing  the  model  are: 

1 .  The  gas  mixture  is  incompressible,  ideal  fluid. 

2.  The  flow  in  fuel  cell  is  laminar  everywhere. 

3.  The  effective  diffusion  coefficient  of  reactant  gases  in  the 
porous  media  is  constant. 

4.  Isothermal  conditions. 

5.  The  diffusion  layer,  catalyst  layer  and  membrane  are 
isotropic  and  homogeneous,  characterized  by  effective 
permeability  and  uniform  porosity. 

6.  Butler- Volmer  kinetics  governs  the  electrochemical  re¬ 
action. 

7.  The  water  produced  by  electrochemical  reaction  at  the 
cathode  side  is  in  vapor  state. 

8.  The  membrane  is  impervious  for  reactants  gases. 

The  governing  equations  are  identical  for  both  conven¬ 
tional  and  interdigitated  flow  fields,  which  will  be  described 
in  the  following  section.  Assuming  that  the  fuel  cell  is  re¬ 
peated  periodically  along  y-direction,  in  order  to  save  the 
computation  source  and  computation  time,  we  take  one  cur¬ 
rent  collector,  half  of  two  flow  channels  around  which  as 
computation  domain,  illustrated  in  Fig.  1. 
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Fig.  1.  Schematic  of  the  flow  fields  structure  and  the  computation  domain  for  a  conventional  flow  field  (a)  and  interdigitated  flow  field  (b).  (c)  A  typical 
cross-section  of  the  domain  and  the  locations  of  different  layers. 


2.1.  Continuity  equation 

The  conservation  of  mass  equation  for  the  gas  mixture 
Vfoow)  =  Sm 


where  Sk  stands  for  the  source/sink  term  of  the  k’th  species 
induced  by  the  electrochemical  reaction  occurred  in  the  ac¬ 
tive  catalyst  layers. 


2.2.  Momentum  conservation 


where  e  is  the  porosity  of  the  porous  media,  which  is  equal 
to  unit  for  the  gas  channels,  p  the  density,  and  u  the  intrin¬ 
sic  fluid  velocity  vector,  su  reflects  the  superficial  velocity 
in  the  porous  media.  Sm  denotes  the  mass  source  term  con¬ 
sidering  the  electrochemical  aspects  of  the  fuel  cell,  thus 
they  are  zero  in  most  of  the  computation  domain  except  for 
the  activated  reaction  sites  as  the  anode  and  cathode  cata¬ 
lyst  layers.  These  terms  corresponding  to  the  consumption 
of  hydrogen  and  oxygen  in  the  anode  and  the  cathode,  and 
the  production  of  water  in  the  cathode  can  be  determined  as 
following: 

Sm  =  J2  V  (2) 

k 


The  fluid  flow  in  the  fuel  cell  can  be  described  by  the 
general  equation  as: 

V(spuu)  =  —  sV  p  +  V(£/xeffVw)  +  Su  (3) 

where  p  denotes  the  pressure,  /xeff  the  effective  viscous  co¬ 
efficient.  Because  the  fluid  flowing  in  the  channels,  backing 
layers  and  catalyst  layer,  membrane  is  different,  gas  mix¬ 
ture  and  liquid,  so,  /zeff  stand  for  gas  viscous  coefficient 
for  gas  mixture  in  the  channel  and  backing  layer,  and  liq¬ 
uid  viscous  coefficient  for  liquid  in  the  catalyst  layer  and 
membrane.  Furthermore,  mass-weighted  mixing  law  gives 
viscosity  of  the  gaseous  mixture.  The  momentum  source  Su 
is  used  to  describe  the  Darcy’s  drag  force  exerted  by  the 
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fluid  on  the  solid  surface  in  the  porous  electrode,  catalyst 
layer  and  membrane.  For  the  fluid  flow  in  the  gas  channel, 
the  equation  becomes  the  Navier-Stokes  equation.  The  mo¬ 
mentum  conservation  in  the  porous  media  can  be  reduced 
to  the  well-known  Darcy’s  law  because  the  inertia  and  vis¬ 
cous  terms  is  so  small,  and  can  be  neglected.  However,  the 
macroscopic  velocity  changes  rapidly  at  the  interface  be¬ 
tween  the  channel  and  gas  diffusion  layer,  the  macroscopic 
viscous  force  and  inertial  force  must  be  considered  in  the 
momentum  equation.  So,  a  general  equation  with  appropri¬ 
ate  source  terms  [1 1]  is  more  suitable  for  describing  the  flow 
in  the  whole  fuel  cell.  The  momentum  source  terms  for  dif¬ 
ferent  layers  can  be  as  [12]. 


J  a  —  Ja,ref 


CH2 

^442,ref 


Ka 


X 


exp  (^4r>7act)  -  exp 


where  rfact  the  surface  overpotential,  is  defined  as: 

OsLCt  —  ^a,c  Vqq  (9) 


0 


/,eff 
l1  2 
- £11 

k 

/4eff  k<u 

~T  Sm£mcll  +  ~j^-ZfCfFV(pm 

/Cp  /Cp 

ueff  K 

——£mu  +  -f-ZfCfFV(pm 

kp  kp 


gas  channel 
backing  layer 

catalyst  layer 

membrane 

(4) 


Where  cpm  stands  for  the  phase  potential  of  the  electrolyte 
membrane,  £m  is  membrane  porosity  and  £mc  is  the  volume 
fraction  of  the  membrane  phase  in  the  catalyst  layer,  k  is 
permeability  of  the  backing  layer.  kp,  k$  stand  for  the  hy¬ 
draulic  permeability  and  electrokinetic  permeability.  Zf,  Cf 
are  fixed  site  charge  and  fixed  charge  concentration.  This 
equation  represents  a  generalized  Darcy  equation  for  the 
backing  layer,  and  a  generalized  form  of  Schlogl’s  equation 
for  the  catalyst  layer  and  membrane. 


where  cp^c  and  cpm  denote  the  potentials  for  the  carbon  phase 
and  membrane  phase,  respectively,  at  the  catalyst  layer.  Voc 
is  the  reference  open-circuit  potential  of  the  electrode. 

2.4.  Potentials 

The  electrical  potentials  of  the  membrane  phase  cpm  and 
carbon  phase  cpa,  <Pc  are  introduced  following  the  concept  of 
mean  parameters.  These  two  types  of  potentials  govern  the 
motion  of  protons  and  that  of  electrons,  respectively.  Taking 
local  electro-neutrality  into  account,  the  electronic  current 
and  ionic  current  produced  or  consumed  in  the  catalyst  layers 
lead  to  a  voltage  drop  via  Ohm’s  law  according  to: 

i  =  ±oWcp  (10) 

The  governing  equation  for  each  sort  of  current  follows  the 
current  conservation: 

V/  =  j  (11) 


2.3.  Species  conservation  equations 

The  species  conservation  equation  for  the  gas  mixture  is 
V(enCt)  =  V(£>f  VQ)  +  Sk  (5) 

Here,  k  denotes  chemical  species  that  include  hydrogen, 
oxygen,  nitrogen  and  water.  is  the  effective  diffusion 
coefficient.  The  S  depending  on  the  transfer  current  can  be 
described  as  following: 


Thus,  the  three  potentials  obey  the  following  equations: 


V  (crm  V  cpm )  — 


V(oraV<pa)  =  - 

V(<7cV<pc)  =  « 


—7a, 

X\  <  X  <  X2 

0 

X2  <  X  <  V3 

(12) 

[  7c, 

V3  <  X  <  X4 

_  a 

O  ^ 

<  X\ 

<  X  <  X2 

(13) 

P  1 

X3  <  X  <  X4 

X4  <  X 

(14) 

A 

IF 


for  H2 


A 

4  F 


for  O2 


A 

IF 


for  H20 


After  obtaining  the  phase  potentials,  the  current  density  in 
the  external  circuit  I  can  be  calculated  by  the  following 
(6)  expression: 


where  j  stands  for  the  transfer  current  density,  which  can  be 
calculated  by  the  following  Butler- Volmer  kinetics  expres¬ 
sion: 


where  H  is  the  height  of  the  catalyst  layer,  L  the  length.  As 
<Pa  Ix=jc0  is  taken  to  be  zero,  then  cpc  |X=X5  -  <pa \X=XQ  =  cpc  \x=x5 
gives  total  voltage  of  the  whole  fuel  cell. 
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2.5.  Boundary  conditions 

The  computational  domain  covers  a  part  of  the  cell  assem¬ 
bly,  which  is  repeated  periodically  along  y-direction.  There¬ 
fore,  on  the  bottom  and  top  ends  (y  =  0  and  y  =  H)  periodic 
boundary  conditions  for  all  variables  are  imposed:  f\y=o  = 
f\y=H ,  where /stands  for  all  potentials,  concentrations,  and 
velocity. 

Boundary  condition  corresponding  to  membrane  phase 
potential  equation  is  zero  flux  along  all  boundaries  of  the 
computational  domain,  giving  Von  Neumann  boundary  con¬ 
dition  dcpm/dn  =  0.  In  the  current  collector  the  carbon  phase 
potentials  are  fixed  (<^  =  0,  (p®).  Solution  of  the  problem  for 
a  given  provides  the  mean  current  density  in  the  fuel  cell, 
i.e.  a  point  on  the  I-V  curve.  The  electron  current  does  not 
flow  to  the  gas  channel  and  the  membrane,  the  normal  com¬ 
ponent  of  electron  current  at  the  channel  surfaces  is  zero: 
dcpa,c/dx  =  0.  And  the  same  condition  is  imposed  on  both 
sides  of  the  membrane. 

As  the  membrane/catalyst  interface  is  assumed  to  be  im¬ 
permeable  to  gases,  the  flux  of  gases  is  zero  at  this  interface 
dCk/dx  =  0.  The  same  condition  is  imposed  on  surfaces  of 
metal  electrodes. 

At  the  inlets,  the  fluid  is  supposed  to  flow  into  the  chan¬ 
nel  at  a  known  velocity,  and  the  species  concentrations  are 
known.  Thus,  the  following  Dirac  boundary  conditions  are 
prescribed: 

U  =  Uin,  V  =  W  =  0,  Ck  =  C*.  in 

At  the  outlets,  the  fully  developed  conditions  are  applied  to 
velocity  and  species  concentration  fields. 

2.6.  Numerical  procedures 

The  conservation  equations,  including  mass,  momentum, 
species  and  charges,  were  converted  to  a  finite-difference 
form  by  the  control  volume  method.  An  orthogonal  non- 
uniform  grid  for  computational  discretion  was  employed  in 
this  work.  The  flow  channels  of  anode  and  cathode  are  di¬ 
vided  into  12x20x360  equally  sized  grid  cells.  The  y-z 
planes  of  backing  layer,  catalyst  layer  and  membrane  are  di¬ 
vided  into  40x360  grids,  however,  in  the  v-direction  through 
the  cell,  backing  layers  are  equally  divided  into  four  parts, 
catalyst  2  and  membrane  4.  The  solutions  of  the  velocity 
component  equations  are  obtained  in  a  staggered  control 
volume.  The  SIMPLEC  (semi-implicit  method  for  pressure- 
linked  equations  consistent)  scheme  was  adopted  for  the 
pressure  correction.  The  species  equations  are  calculated  fol¬ 
lowing  the  velocity  vector,  coupled  with  calculation  of  the 
potentials  by  electrochemical  reaction  kinetics.  The  unified 
domain  method  was  adopted  for  the  coding  in  this  work;  we 
use  proper  boundary  for  the  parameters,  rather  than  bound¬ 
ary  conditions  for  the  dependent  variables  at  the  interface 
between  different  layers  [13].  The  coupled  set  of  equations 
was  solved  iteratively  until  the  relative  error  in  each  field 
achieved  the  specific  convergent  standard  (usually  10-5). 


3.  Results  and  discussions 

The  results  presented  below  were  obtained  for  the  param¬ 
eters  and  base  conditions  listed  in  Table  1 ,  which  is  obtained 
mainly  from  Bernardi  et  al.  [14]. 

Fig.  2  illustrates  two-dimensional  velocity  vectors  in  the 
cathode  x-y  plane  at  z  =  L/2  for  both  conventional  and 
interdigitated  flow  fields.  It  can  be  seen  that  the  magnitude 
of  velocity  in  the  interdigitated  flow  field  is  one  or  two  or¬ 
ders  larger  than  that  in  the  conventional  flow  field.  Fig.  2a 
shows  the  velocity  vectors  in  the  conventional  flow  field, 
these  velocities  are  pointing  downward  from  the  membrane 
to  the  cathode  channel.  It  is  because  that  in  the  cathode 
oxygen  consumed  on  the  membrane  surface  is  less  than  the 
water  vapor  produced  on  the  membrane  surface.  Fig.  2b 
displays  the  velocity  field  in  the  interdigitated  flow  field. 
It  is  evident  that  there  exist  bypass  flow  around  the  cur¬ 
rent  collector  shoulder  and  also  strong  forced  convection 
along  the  reaction  surface.  The  flow  from  the  inlet  chan¬ 
nel  to  the  outlet  channel  can  enhance  providing  the  reac¬ 
tants  to  the  catalyst  layers  and  removing  the  products  from 
the  catalyst  layers.  Thus,  the  performance  of  the  fuel  cell 
can  be  improved  essentially,  especially  at  higher  current 
density. 

Table  1 


Physical  parameters  and  base  conditions 


Quantity 

Value 

Operating  temperature  (°C) 

80 

Pressure  of  anode  channel  (atm) 

1 

Pressure  of  cathode  channel  (atm) 

1.5 

Length  of  the  channel,  L  (cm) 

6.0 

Channel  height  (cm) 

0.05 

Fuel  cell  height  (cm) 

0.2 

Channel  width  (cm) 

0.07 

Backing  layer  width  (cm) 

0.024 

Catalyst  layer  width  (cm) 

0.006 

Membrane  width  (cm) 

0.025 

«a 

0.5 

«c 

1.0  [16] 

Ka 

0.5 

Yc 

1.0 

Hydraulic  permeability  of  membrane,  kp  (cm2) 

i.8  x  icr14 

Permeability  of  backing  layer,  k  (cm2) 

1.76  x  10“7 

Electrokinetic  permeability,  k ^  (cm2) 

7.18  x  10“16 

Hydrogen  mole  fraction  at  the  inlet  of  the  cathode 

1.0 

Membrane  porosity,  sm 

0.5 

Volume  fraction  of  membrane  in  the  catalyst  layer, 

0.5 

fcmc 

Porosity  of  the  backing  layer,  £o 

0.5 

Membrane  phase  conductivity,  crm  (£2_1  cm-1) 

0.034 

Carbon  phase  conductivity,  crac  cm-1) 

40 

Fixed  charge  concentration,  Cf  (mol cm-3) 

1.2  x  1CT3 

Oxygen  mole  fraction  at  the  inlet  of  the  cathode 

0.21 

Nitrogen  mole  fraction  at  the  inlet  of  the  cathode 

0.79 

Reference  transfer  current  density  at  anode,  ya,ref 

5.0  x  102 

(A  cm-3) 

Reference  transfer  current  density  at  cathode, jcjef 

4.4  x  10“4 

(A  cm-3) 

Velocity  at  the  inlet  (cms-1) 

50 
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0.0002  m/s 


0.02  m/s 


Fig.  2.  Two-dimensional  velocity  fields  in  the  cathode  x-y  plane 
cross-backing  layer  and  flow  channels  at  z  =  L/2  for  V^n  =  0.7  V:  (a) 
conventional  flow  field;  (b)  interdigitated  flow  field. 

According  to  species  transport  equation  5,  the  species  flux 
is  supported  by  two  processes:  convection  and  diffusion. 
Fig.  3  shows  the  ratio  of  both  absolute  components  of  the 
flux  of  the  oxygen  in  the  cathode  x-y  plane  at  z  =  L/2  for 
the  conventional  and  interdigitated  flow  fields.  It  is  seen  that 
in  the  conventional  flow  field  the  convection  term,  \suC\,  is 
two  or  three  order  less  than  the  diffusion  term,  |DVC|,  the 
reactants  is  mainly  transported  by  the  diffusion  mechanism, 
as  displayed  in  Fig.  3a.  This  is  because  that  the  velocity  for 
the  conventional  flow  field  is  very  small,  however,  the  gra¬ 
dient  of  concentration  is  higher.  However,  for  the  interdigi¬ 
tated  flow  field,  the  forced  flow  enhance  convection  flux  in 
the  diffusion  layer,  the  convection  term  take  more  important 
effect  to  the  transport  process,  and  enhance  the  gas  mixture 
transport,  as  shown  in  Fig.  3b.  It  also  indicates  those  inter- 


(a)  at  z=L/2 


(b)  at  z=L/2 

Fig.  3.  Ratios  folfc  of  absolute  values  of  diffusion  /b  =  \DWC\  and 
convective  fc  =  \suC\ ,  components  of  oxygen  flux  in  the  cathode  x-y 
plane  at  z  =  L/2  for  Vcen  =  0.7V:  (a)  conventional  flow  held;  (b) 
interdigitated  how  held. 

digitated  flow  fields  effectively  convert  the  transport  of  the 
reactant  or  product  gases  to  or  from  the  catalyst  layers  from 
a  diffusion  mechanism  to  a  forced  convection  mechanism. 

Fig.  4  shows  the  water  vapor  mole  fraction  in  the  y-z 
plane  at  the  interface  between  the  cathode  catalyst  layer  and 
the  porous  backing  layer  for  the  conventional  and  interdig¬ 
itated  flow  fields.  From  Fig.  4a,  in  the  conventional  flow 
field,  the  water  vapor  distribution  is  symmetric  around  the 
centerline  in  z-direction  for  the  same  flow  in  the  channel, 
and  increase  along  the  channel  for  the  electrochemical  re¬ 
action.  In  the  interdigitated  flow  field,  however,  maximum 
water  vapor  concentration  appears  in  the  outlet  portion  of 
the  outlet  channel,  because  the  flow  from  inlet  channel  into 
electrode,  and  out  from  the  outlet  channel,  as  illustrated  in 
Fig.  4b. 

Fig.  5  displays  the  local  current  density  distribution  in  the 
cathode  catalyst  layer  (in  the  y-z  plane)  with  voltage  0.7  V 
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Fig.  4.  Water  mole  fraction  distributions  at  the  interface  between  the  catalyst  layer  and  the  membrane  for  Vcen  =  0.7  V:  (a)  conventional  flow  field;  (b) 
interdigitated  flow  field. 


for  the  conventional  and  interdigitated  flow  fields.  The  lo¬ 
cal  current  density  has  similar  trends  because  it  depends  on 
the  oxygen  amount  transported  to  the  catalyst  layer.  Fig.  5a 
shows  the  local  current  density  distribution  in  the  conven¬ 
tional  flow  field.  The  current  density  distribution,  which  is 
controlled  by  oxygen  diffusion,  is  symmetric  around  the 
centerline  in  z-direction  for  the  same  flow  in  the  channel. 


The  local  current  density  in  the  front  of  the  current  collec¬ 
tor  is  smaller  than  that  in  front  of  the  channel,  because  the 
limitation  of  the  diffusion,  the  reaction  surface  cannot  be 
applied  fully.  The  current  density  in  the  interdigitated  flow 
field  is  controlled  by  forced  convection,  shown  in  Fig.  5b. 
Near  the  entrance  of  the  flow  channel  in  the  interdigitated 
flow  field,  the  local  current  density  is  higher  because  the 
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Fig.  5.  Local  current  density  distributions  in  the  y-z  plane  for  Vceii  =  0.7  V: 
(a)  conventional  flow  field;  (b)  interdigitated  flow  field. 


oxygen  concentration  gradients  across  the  backing  layer  is 
larger,  followed  by  diminishing  local  current  density  due  to 
the  reduced  oxygen  concentration  along  the  flow  channel 
direction  for  the  electrochemical  reaction. 

Fig.  6  shows  the  pressure  drop  for  the  conventional  and 
interdigitated  flow  fields.  The  interdigitated  flow  field  causes 
larger  pressure  drop  than  the  conventional  flow  field  for  its 
fluid  flowing  through  the  porous  electrode. 

The  comparison  of  performance  for  conventional  and  in¬ 
terdigitated  flow  fields  is  illustrated  in  Fig.  7.  It  is  obvious 
that  the  interdigitated  flow  field  outperforms  the  conven¬ 
tional  flow  field,  especially  at  high  current  densities.  The 
limit  current  density  for  a  fuel  cell  with  the  interdigitated 
flow  field  is  improved  one  time  corresponding  to  the  fuel 
cell  with  the  conventional  flow  field.  This  is  because  that 
high  oxygen  transfer  rate  is  required  to  be  provided  at  high 
current  density,  which  cannot  be  satisfied  by  the  conven¬ 
tional  flow  field  for  the  diffusion  limitations,  whereas  can 


Fig.  6.  Pressure  drops  in  the  conventional  and  interdigitated  flow  fields 
for  different  current  density. 


Fig.  7.  Comparison  of  performance  of  fuel  cell  with  conventional  and 
interdigitated  flow  fields. 

be  fulfilled  by  the  interdigitated  flow  field  by  strong  forced 
convection. 

The  availability  of  our  model  and  its  simulation  was  ver¬ 
ified  by  comparing  estimated  and  experimental  polarization 
curves  obtained  by  Nguyen  [15],  as  displayed  in  Fig.  8.  The 
numerical  model  has  been  run  with  parameters  reported  in 
Table  1  and  the  conditions  similar  to  experimental  case,  i.e. 
temperature  is  298.15  K,  and  pressure  is  1  atm.  It  indicates 
that  the  predicted  results  agree  well  with  experimental  data. 
It  should  be  pointed  out  that  the  reference  transfer  current 
density  was  adjusted  to  achieve  the  good  agreement  between 
the  present  model  predictions  and  the  experimental  data. 
Both  the  modeling  results  and  the  experimental  data  show 
the  same  factor  of  the  fuel  cell  performance  increase  when 
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Fig.  8.  Comparison  between  experimental  (scattered  dots)  and  estimated 
polarization  curves  (full  lines):  (■)  interdigitated  flow  field;  (T)  conven¬ 
tional  flow  field. 

interdigitated  flow  field  is  used.  Especially,  an  interdigitated 
flow  field  can  greatly  enhance  limiting  current  density,  be¬ 
cause  mass- transport  limitations  become  main  hindrance  for 
performance  under  high  current  density. 

4.  Conclusions 

A  steady-state,  three-dimensional  mathematical  model 
was  employed  to  investigate  and  compare  the  characteristics 
of  a  proton  exchange  membrane  fuel  cell  with  the  interdig¬ 
itated  and  conventional  flow  fields.  The  multidimensional 
flow  structure,  distributions  of  reactants  concentrations  and 
current  densities  for  the  fuel  cell  has  been  calculated.  Com¬ 
parison  with  a  conventional  flow  fields  design  fuel  cell,  an 
interdigitated  flow  fields  can  provide  high  reactants  to  cat¬ 
alyst  and  products  from  catalyst;  enhance  limiting  current 
density  and  polarization  performance  by  forced  convec- 
tional  flow  fields  through  the  electrode.  However,  there 
have  larger  pressure  drop  between  the  inlet  and  the  outlets 
in  the  fuel  cell  with  interdigitated  flow  field,  which  may 
be  a  main  reason  hindering  interdigitated  design  applica¬ 
tion.  Polarization  curves  predicted  by  numerical  modeling 
compared  well  with  experimental  data  in  the  published 
literature. 
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